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The zero-pressure zero-temperature equilibrium volumes and bulk moduli are calculated for the 
light actinides Th through Pu using two independent all-electron, full-potential, electronic-structure 
methods: the full-potential linear augmented-plane-wave (FLAPW) method and the linear combi- 
nations of Gaussian-type orbitals-fitting function (LCGTO-FF) method. The results produced by 
00 ' these two distinctly different electronic-structure techniques are in excellent agreement with each 

other, but differ significantly from previously published calculations using the full-potential linear 
muffm-tin-orbital (FP-LMTO) method. The theoretically calculated equilibrium volumes are in 
some cases nearly 10% larger than the previous FP-LMTO calculations, bringing them much closer 
to the experimentally observed volumes. We also discuss the anomalous upturn in equilibrium 
volume seen experimentally for a-Pu. 
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I. INTRODUCTION 

The light-actinide metals, from Th to Pu, pose a severe challenge to modern electronic-structure theory due in part 
to the existence of highly directional and high density-of-states /-electron bonding, which promotes the formation of a 
large number of exceptionally complicated crystal structures. In fact, the light actinide metals (excluding Ac) together 
with Ce are the only elemental solids known to exhibit /-bonding under ambient conditions. In addition, the heavy 
nuclei of the light-actinide metals induce large relativistic effects on their valence-band structures, and usually require 
inclusion of spin-orbit-coupling effects. Despite these difficulties, many all-electron electronic-structure calculations 



have been carried out for the light-actinide metals.ETO To the best of our knowledge, only ane full-potential technique, 
the full-potential linear muffin-tin-orbital (FP-LMTO) method has thus far been applied,cl but taken all together, the 
< , cumulative theoretical results for the bulk properties (such as the equilibrium volume and bulk modulus) are at best 
7-H ' ambiguous, even for the simplest of these materials, fee Th. 

Local-density approximation (LDA)I3 calculations of the ambient properties of fee Th using the linear muffin-tin- 
orbital (LMTO) method within the atomic-sphere approximation (ASA) produced zero-pressure volumes that range 
from slightly compressed relative to experiment to, slightly expanded, depending on the inclusion or neglect of spin- 

S orbit corrections and_combined-correction terms .EJ El When those calculations were repeated using a full-potential 
(FP-LMTO) method,Q the presumably more accurate FP-LMTO method yielded a volume nearly 20% smaller than 
, the best LMTO- AS A result and experiment. If this FP-LMTO result is correct, fee Th exhibits the largest LDA- 
induced contraction for any elemental solid other than fee Pu (and the other heavy fermion systems), which will be 
discussed later, and the largest volume discrepancy between FP-LMTO and LMTO-ASAjpalculations found to date. 
O ■ Additional FP-LMTO calculations using the generalized-gradient-approximation (GGA)El produced a zero-pressure 
volume 10% larger than the FP-LMTO LDA result, but still nearly 10% smaller than the measured value, which is 
another rather anomalous result. In fact, FP-LMTO GGA calculations have consistently yielded volumes that are 
significantly smaller than experiment for.-all of the light-actinide metals, with the size of the error ranging from just 
under 6% for Np up to about 10% for Th.Eil Recent full-charge-density (FCD)-LMTO calculations on Th found an LDA 
lattice constant that lies between the LDA values found in the earlier LMTO-ASA and FP-LMTO calculations and a 
GGA lattice constant that is significantly larger than the FP-LMTO GGA result. Taken together, these inconsistent 
results raise serious doubts about the reliability of the existing theoretical calculations for the light-actinide metals. 

To resolve these issues, we have carried out electronic-structure calculations on the light actinides using two, in- 
dependent methods: the linear combinations of Gaussian-type orbitals-fitting function (LCGTO-FF) methodj3 as 
implemented in the program GTOFFEl and the full-potential linearized augmented plane-wave (FLAPW) method.EJ 
Although both methods are all-electron, full-potential techniques capable of producing high-precision total energies 
within the LDA and GGA approximations, they employ radically different numerical approximations and basis sets. 
Application of both methods to a single system thus provides assurance that the final results obtained are free of 
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numerical artifacts. For this investigation, LDA and GGA calculations were carried out scalar-relativistically and 
fully-relativistically (i.e., with spin-orbit effects included). The FLAPW method includes spin-orbit coupling as a 
perturbation,!!!! while the LjGGTO-FF performs a non-perturbative treatment of spin-orbit effects, which has only 
recently been implcmcnted.tj It is of particular interest to compare the scalar and fully relativistic results across 
the light actinides; not only does this comparison justify the perturbative approach to spin-orbit coupling, it also 
illustrates the effects of spin-orbit coupling on the bulk properties. 

Section JfL-briefly recapitulates the FLAPW and LCGTO-FF methods, which have been discussed extensively 
elsewhere, BrL3 along with the details of the present calculations. Section |l| then discusses the results for equilibrium 
volumes and bulk moduli for Th through Pu, both in the fee (Th— >Pu) and a (Pa^Pu) crystal phases (when different 
from fee). We have included the fee phase for Pa^Pu as it has often been used as a surrogate for the much more 
complex a structures. Section IV presents our concluding remarks. 



II. METHODS 



The LCGTO-FF technique is distinguished from other electronic structure methods by its use of three independent 
GTO basis sets to expand .the orbitals, charge density, and exchange-correlation, (XC) integral kernels; here using 
either the Hedin-LundqvistB LDA model or the Perdew-Wang 91 GGA model.0 The charge-fitting functions are 
used to reduce the total number of Coulomb integrals required by replacing the usual 4-center integrals in the total 
energy and one-electron equations with 3-center integrals. The_jcharge-fitting function coefficients are determined 
by minimizing the error in the Coulomb energy due to the fit;E£l thereby allowing high-precision calculations with 
relatively small basis sets. The least squares XC fit used here acts as a sophisticated numerical quadrature scheme 
capable of producing accurate results on a rather coarse numerical integration mesh. Scalar-relativity.^ras recently 
implemented and tested p,GTOFFll3 using the nuclear-only, Douglas-Kroll-Hess (DKH) transformatiorotJ developed 
by Habeclen and R6sch,ll3 and scalar relativistic cross-product terms and spin-orbit effects have also been included 
recently]!!! thus enabling a fully relativistic treatment. 

The precision of any LCGTO-FF calculation will, of course, be largely determined by the selection of the three GTO 
basis sets. In this work, 23s20pl5dll/ uncontraeted orbital-basis sets were derived for the light actinides from atomic 
basis sets tabulated by Minami and MatsuokaO These basis sets were reduced to 17sl4pll(i7/ contracted basis sets 
using coefficients from atomic calculations that used the same density functional model as the crystal calculations. 
The charge density and XC integral kernels were fitted with 25s and 21s basis sets, respectively, selected on the basis 
of prior experience with LCGTO-FF calculations. All of these basis sets can be obtained from the authors. 

All Brillouin zone (BZ) integrations required for the LCGTO-FF calculations were carried out on a uniform mesh 
with 72 irreducible k-points in the fee Brillouin zone (163 for a-Pa and 63 for a-U), using a Gaussian-broadened, 
histogram, integration technique, with a broadening factor ranging from 4 mRy for Np to 10 mRy for Th. The SCF 
cycle was iterated until the total energy varied by less than 0.02 mRy/atom. _ . 

The all-electron FLAPW method that we have used has been described by Singh;L3 here we discuss only the essential 
features relevant for comparison the with other full-potential techniques. Energy parameters for the FLAPW basis 
functions were set near the centers of their respective bands by monitoring the eigenvalues of the calculation whose 
volume lay closest to equilibrium. Local orbitala!3 were added to enhance the variational freedom and allow the 
semi-core 5d, 6s, and 6p orbitals to be treated along with the valence electrons, with the added energy parameter 
used to simultaneously treat the residual s and p character of the valence electrons. Local orbitals with / character 
were also used to maximize valence-band accuracy and ensure orthogonality with the core 4/ states. The results 
that we obtained were found to be insensitive to small changes in all energy parameters. Spin-orbit coupling was 
incorporated for the valence electrons at the second- variational level,t2l in which the effects of spin-orbit coupling are 
treated perturbatively in a set of scalar-relativistic orbitals found within an energy window of specified width. The 
results obtained here were found to be stable with respect to the size of this window. 

The size of the FLAPW basis is determined by a plane- wave cutoff, K max , whose value was given by the relation 
RMTK max = 9.0, which was found to be satisfactory in determining the bulk properties. For BZ integrations, 60 
irreducible points in the fee Brillouin zone (80 poiats in a-Pa, 126 in a-U, 8 in a-Np, and 16 in a-Pu) were sampled 
according to the scheme of Monkhorst and Packjl3 and all BZ integrations were reduced to simple sums of a Fermi- 
Dirac-function temperature-broadened integrand, with an effective temperature of 2.0 mRy. Increasing the number 
of points sampled in the BZ changed the total energies by less than 0.25 mRy/atom for the a phases, and 0.75 
mRy/atom for the fee structures. The values thus obtained for the equilibrium volumes and bulk moduli were found 
to be converged with respect to the number of BZ points sampled. All self-consistent calculations were iterated 
until the total energy changed by less than 0.01 mRy/atom. The same LDA and GGA models were used as in the 
LCGTO-FF calculations. 
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As is evident from the preceding discussion, the FLAPW and LCGTO-FF methods differ significantly. First, the 
FLAPW method utilizes a diffuse orbital-basis set that can be systematically enriched by simply adding more plane 
waves. This is in sharp contrast to the local basis set used in the LCGTO-FF method, which tends to be somewhat 
more difficult to improve systematically. On the other hand, the FLAPW basis functions have discontinuities in their 
second derivatives at the muffin-tin sphere boundaries that do not exist in the LCGTO-FF basis sets. 

The implementations of relativity used by the LCGTO-FF and FLAPW methods also differ dramatically. The 
FLAPW method avoids the well-known variational collapse problem by using a basis set that is formed from pure 
electronic solutions to the muffin-tin Dirac-Kohn-Sham equations inside the muffin-tin spheres. (In the interstitial 
region, non-relativistic solutions are used.) The non- muffin-tin corrections to the potential are then accounted for 
perturbatively at each iteration. This approach to solving the full-potential Dirac-Kohn-Sham equations should be 
accurate provided the muffin-tin solutions do not differ greatly from the full-potential solutions and the effects of 
relativity are small in the interstitial region. On the other hand, the LCGTO-FF method decouples the electron 
and positron degrees of freedom of the full-potential Dirac-Kohn-Sham equations to second-order in the ratio of the 
effective potential to the combined kinetic and rest mass energies. In general, this approach will not be as accurate 
near the nucleus as the FLAPW method, but will be more accurate in the interstitial region. While spin-orbit effects 
are tcatcd perturbatively in the FLAPW method, the LCGTO-FF method doubles its full Hamiltonian and overlap 
matrices. A comparison between LCGTO-FF and FLAPW scalar and fully relativistic results thus provides a check 
on the second-variational treatment used in the FLAPW method for the light actinides. 

We wish to emphasize that that the LCGTO-FF and FLAPW techniques may be viewed as complementary methods, 
to the extent that they utilize different numerical approximations. Thus, if one of the two methods were to encounter 
difficulties due to inadequacies in its basis set or implementation of relativity, it is unlikely that the second method 
would exhibit the same problem. This makes a joint study of the type carried out here especially valuable for 
investigating systems, like the light actinides, that have previously exhibited a significant sensitivity to computational 
details. 



III. RESULTS 

For each choice of system (Th^Pu), method (LCGTO-FF or FLAPW), model (LDA or GGA), and level of 
relativity (scalar or full), total energies were calculated for five or six volumes lying near the energy minimum. The 
calculated energies for each combination were then fitted with either a cubic function of the volume or a second- 
order Birch fit $3 to obtain the zero-pressurej-wplupre (Vq) and bulk modulus (-Bo). The fitted results are compared 
with previous calculationstTH and experimentnJ~Ed in Tables || and [[J. We have performed these calculations in both 
the fee (Pearson cF4) and experimentally determined a structures, in order to determine to what extent the cubic 
close-packed structure can serve as a surrogate for the more complex (and much more computationally demanding) a 
structures. No attempt was made to relax the structural parameters of the a phases in this study. Indeed, the a-Np 
and a-Pu phases were already too demanding for the current implementation of the LCGTO-FF method. 

Inspection of the results for the equilibrium volumes listed in Table | reveals several interesting features. First, 
the LCGTO-FF and FLAPW methods give nearly identical scalar and fully relativistic results for Th— >Pu. This 
level of agreement between such disparate methods provides a high degree of confidence in the quality of the results 
produced by both methods. In contrast, the fully-relativistic FLAPW and FP-LMTO resultso differ substantially for 
both the LDA and GGA models, with the volumes differing by roughly 3-10% in each case, with the notable exception 
of a— Pu. This is an exceptionally large disagreement for two full-potential methods. The good agreement between 
the LCGTO-FF and FLAPW results lends confidence to the FLAPW calculations. In addition, the fully-relativistic 
FP-LMTO GGA volume for Th is anomalously small compared to experiment. Taken together, these features strongly 
suggest that the sourcer,of the disagreement between the present results and the FP-LMTO results is some error in 
the earlier calculations. □ 

The FLAPW and LCGTO-FF GGA results are compared with experiment (in most cases only room temperature 
experimental data is available) in Figured for the a phases. The agreement between theory and experiment is quite 
good. A similar plot is shown in Figure |ffor the fee crystal structures, in which we see that the fee have significant 
errors when compared to the a phase experimental volumes. The fee phases do, however, obey the same general 
trend with regard to spin-orbit contributions; spin-orbit effects have negligible impact on the equilibrium volumes for 
Th^U, and cause an expansion of the lattice for Np and Pu. Similarly, GGA results are always expanded relative to 
LDA. 

The upturn in equilibrium volume for Pu seen in Figure [l] has been a subject of some contention. Based on 
LMTO-ASA studies, this upturn was originally ascribed to spin-orbit effects on the valence electrons .Em FP-LMTO 
studies by Wills and Eriksson,Q however, not only did not find an upturn, but found a contraction with the inclusion 



3 



LA-UR-99-2778 



February 7, 2008 



submitted to Physical Review B 



of spin-orbit coupling, leading them to hypothesize that the LMTO-ASA calculations reached the wrong conclusion 
about the-cffects of spin-orbit coupling due to their treating the 6p semi-core states as core states. Recently, however, 
PenicaudE3 has performed a fully relativistic LMTO-ASA calculation for the light actinides in which the 6p states 
were treated as valence states (in a relativistic j, K-basis), and found no such contraction. Indeed, the LMTO-ASA, 
studies show a systematic trend towards an expansion in volumes when spin-orbit effects are treated perturbativelyu 
and non-perturbatively.E3 Our present studies are consistent with the LMTO-ASA results, and Figure [l] shows that 
spin-orbit effects on the valence electrons are responsible for a gradual increase of the equilibrium volume as one 
proceeds along the sequence of the light actinides, but it does not appear to account for the anomalous upturn in 
volume for a-Pu. 

Another possibility for the upturn is that it represents a finite-temperature effect (since experimental volumes are 
measured at room temperature). Whilq-4|t is true that thermal expansion increases for Np and Pu, however, if the 
measured thermal expansion coefficients^ (resulting in an increase for the volume of 2.4% for Np, and 4.8% for Pu) 
are used to correct the zero-temperature calculations, the two atomic volumes will be the same at room temperature 
(about 130 in atomic units) and not have the expected upturn. It is quite interesting that our density-functional 
calculations predict the atomic volumes so well until we reach Pu. For fee (5-phase Pu, the theoretical volume is too 
small by about 20%. Although this is an exceptionally large error for a GGA calculation, it has been known for many 
years that S-Pu is anomalous due to its position on the boundary between the light actinides that have itinerant 
5/-electrons and the heavy actinides that have localized 5/-electrons (which form a second rare-earth like series). 
Although the volume found here for 5-Pvl is 10% larger than that found with FP-LMTO methodjj that increase is not 
large enough to remove the discrepancy between theory and experiment. However, the large experimental difference 
in volume between the low-temperature a and much higher temperature S phases suggest that new physics (beyond 
LDA or GGA) is responsible for the large expansion of the 5 phase (the enormous difference cannot be explained by 
any known thermal expansion mechanism). It is likely that large intra-atomic electron-electron coulomb correlations 
(a large effective Hubbard U) are responsible for the anomalous properties of <5-Pu. It may even be possible that some 
localization is already apparent in the a phase for Pu, thus accounting for the slight (around 5%) discrepancy between, 
the computed and measured atomic volume, and causing the upturn. In heavy fermion systems, it has been arguecEj 
that the effects of strong electron-electron correlations (which is responsible for the huge specific-heat enhancements) 
leads to heavier band masses (renormalized bands with a reduced effective band width) . In terms of a Friedel model 
description of bondingp3 such an effect (a narrower / band width) would reduce the strength of the / bonding, which 
should expand the lattice. 

Table |fj lists our results for the bulk modulus for the light actinides Th— >Pu, in both the fee and a phases. 
Since the FP-LMTO calculations found a substantial decrease in the equilibrium volumes, it not very surprising that 
they predict a correspondingly larger bulk modulus. Comparison with the experimentally determined bulk moduli, 
however, which is shown in Figure |[ show that our most accurate fully relativistic GGA FLAPW results are in 
somewhat strong disagreement, particularly for the heavier actinides, Np and Pu. We do not recover the turnover 
seen experimentally past Pa. There are several possible reasons for this discrepancy. One is that the bulk modulus 
is more sensitive than the equilibrium volume to numerical errors in our calculations. Another explanation is that 
temperature effects (all of the available experimental data is at room temperature or above) may play an increasingly 
important role, softening the bulk modulus, especially for Np-and Pu. Strong anharmonic effects have been found 
recently in the Debye- Waller temperatures for U, Np, and PuO 

IV. CONCLUSIONS 

In summary, it has been demonstrated that the LCGTO-FF and FLAPW electronic-structure methods produce a 
zero-pressure volume and bulk modulus for the light actinides Th^Pu that are in good agreement with each other but 
differ significantly from previous full-potential results obtained from the FP-LMTO methodjj with the present volumes 
being roughly 3-10% larger. Given the good agreement between the LCGTO-FF and FLAPW methods, it seems likely 
that the earlier FP-LMTO results were in error. If so, this might account for the 3% to 10% underestimate of the 
volumes found for all of the light-actinide crystals with the FP-LMTO method using the GGA exchange-correlation 
potentials. The present results do not affect the long standing problem of the anomalous volume of <5-Pu, which 
presumably is due to correlation effects that are not adequately accounted for in either the LDA or the GGA. 
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TABLE I. The equilibrium volume (atomic units) obtained here for the light actinides with LCGTO-FF and FLAPW 
calculations using the LDA and GGA models, with and without spin-orbit (SO) effects included, are compared to results from 
LMTO-ASA-icalculations with combined correction terms and calculated temperature expansion corrections,!!! FCD-LMTO 
calculations E and FP-LMTO calculations.^! 



Method 



LDA 



cF4 (fee) 
-SO GGA 



+SO 



LDA 



a-phase 
+SO GGA 



-SO 



cF4 



Th 



LMTO-ASA 232.1 229.4 

FCD-LMTO 212.7 233.7 

FP-LMTO 182.1 

LCGTO-FF 200.5 200.4 216.9 

FLAPW 199.7 204.1 219.3 

Experiment 221.7(298K) a 

cF4 



199.9 
216.2 
218.1 



Pa 



tI2 



LMTO-ASA 

FCD-LMTO 

FP-LMTO 

LCGTO-FF 

FLAPW 

Experiment 



171.7 

159.6 
160.2 



149.9 
160.6 
160.7 



182.8 

171.5 
172.3 



160.2 
172.8 
172.8 



178.8 
164.8 

155.9 
156.8 



U 



cF4 



177.5 

175.8 

147.7 

168.7 

155.8 169.3 
168.30(298K) b 

oC4 



157.6 
170.4 
168.4 



LMTO-ASA 

FCD-LMTO 

FP-LMTO 

LCGTO-FF 

FLAPW 

Experiment 



147.0 

136.1 
136.5 



129.4 
138.8 
136.9 



157.1 

146.1 
147.5 



138.6 
149.7 
148.7 



144.4 
131.4 

128.9 
127.9 



Np 



147.1 
123.7 
128.5 



cF4 



141.2 

138.4 
137.7 
138.89(4.2K) C 

oP8 



131.5 
139.7 
140.1 



LMTO-ASA 

FCD-LMTO 

FP-LMTO 

LCGTO-FF 

FLAPW 

Experiment 



137.4 

121.1 
120.8 



116.7 
125.7 
126.9 



144.1 

131.6 
131.4 



125.8 
138.2 
137.9 



125.5 
117.8 



116.4 
Pu 



128.2 
112.0 
118.3 



cF4 



131.8 



124.6 
129.9(293K) C 

mP16 



122.1, 124.2 a 
127.7 



LMTO-ASA 

FCD-LMTO 

FP-LMTO 

LCGTO-FF 

FLAPW 

Experiment 

Reference [20]. 
b Reference [21]. 
c Reference [22]. 
d Reference [33]. 
"Reference [23]. 
'Reference [24]. 
s Reference [25]. 



114.7 
132.8 

111.5 
111.9 



130.9 

109.2 
119.8 
120.2 



140.0 



121.2 
122.3 



168.0(653K) f 



119.2 
134.4 
133.4 



114.7 
112.2 



109.7 



130.9 



131.4 



114.2 117.2 
134.95(294 K) g 



124.2 d 
124.4 



G 
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TABLE II. The bulk modulus (GPa) obtained here for Th-^Pu with LCGTO-FF and FLAPW calculations using the LDA, 
and GGA models, with and without spin-orbit (SO) effects included, are compared to results from FP-LMTO calculationsu 
and room temperature data. 
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f Reference [30]. 
s Reference [31]. 
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FIG. 1. The equilibrium volume of the light actinides in the a-phases as calculated by FLAPW, both in the scalar-relativistic 
(open squares) and fully-relativistic (solid squares) treatments. These numerical results were done using the GGA form of 
exchange and correlation. Experimental results are also shown (solid circles). 
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FIG. 2. The equilibrium volume of the light actinides in the fee phase as calculated by FLAPW, both in the scalar-relativistic 
(open squares) and fully-relativistic (solid squares) treatments. These numerical results were done using the GGA form of 
exchange and correlation. Experimental results are also shown (solid circles), but for the respective a phases. 
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FIG. 3. The bulk moduli of the light actinides in the a phases, as calculated by FLAPW, both in the scalar-relativistic (open 
squares) and fully-relativistic (solid squares) treatments within the GGA, compared to experiment (solid circles). 
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FIG. 4. The bulk moduli of the light actinides in the fee phase, as calculated by FLAPW, both in the scalar-relativistic 
(open squares) and fully-relativistic (solid squares) treatments within the GGA. For comparison, these values are compared to 
experiment (solid circles) in the respective a phases. 
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